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1.  Introduction 

The  last  decade  has  witnessed  a  great  surge  of  academic  interest 
in  the  dynamic  behavior  and  stability  of  mechanical  systems  with  fol¬ 
lower  (circulatory)  forces,  i.e.  forces  which  are  not  derivable  from  a 
potential,  as  evidenced  by  references  [1]  to  [4].  It  is  a  peculiar 
common  feature  of  much  published  analytical  work  on  this  subject  that 
the  possible  physical  origin  of  such  forces  is  not  mentioned.  The 
follower  forces  are  introduced  into  the  analysis  either  through  a 
sketch,  with  forces  being  merely  indicated  by  arrows,  or  through  a 
specified  functional  dependence  of  the  forces  on  generalized  coordi¬ 
nates.  Thus  the  purportedly  physical  problem  is  reduced  immediately 
to  mathematical  analysis  and  the  relationship  to  mechanics  becomes 
most  tenuous.  The  motivation  for  much  of  this  type  of  work  appears 
to  have  been  sheer  curiosity  in  determining  the  sometimes  unexpected 
behavior  of  an  imagined  system,  rather  than  a  modeling  and  an  explana¬ 
tion  of  observed  phenomena. 

One  of  the  possible  ways  to  realize  follower  forces  is  by  convey¬ 
ing  fluid  through  articulated  or  continuously  flexible  pipes.  The 
pioneering  work  in  this  class  of  mechanical  systems  with  follower 
forces  was  carried  out  by  Benjamin  [5].  He  examined  analytically  and 
experimentally  a  system  with  two  degrees  of  freedom  consisting  of  two 
articulated  pipe  segments,  constrained  to  motion  in  a  plane  like  a 
double  pendulum,  and  conveying  fluid.  It  was  found  that,  depending 
upon  the  values  of  relevant  parameters,  the  system  could  lose  stability 
by  static  buckling  (passage  to  ctn  adjacent  equilibrium  position, 
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divergence)  or  by  dynamic  instability  (oscillations  with  increasing 
amplitude,  flutter).  Several  other  workers  considered  various  aspects 
of  the  dynamic  behavior  of  continuous  pipes  conveying  fluid,  as  men¬ 
tioned  in  [4].  However,  the  transition  from  one  type  of  loss  of  stabil¬ 
ity  to  another  was  not  investigated  systematically  and  quantitatively. 

The  purpose  of  the  present  investigation  is  to  study  the  various 
modes  of  instability  of  a  mechanical  system  whose  non-conservative 
character  arises  in  a  completely  natural  manner,  and  to  verify  that 
the  transitions  between  the  instability  modes  (as  predicted  by  the 
linearized  dynamic  equations  of  motion)  are  indeed  experimentally  ob¬ 
tainable.  To  this  end,  a  spatial  system  of  articulated  pipes  is  con¬ 
sidered.  This  system  is  a  generalisation  of  the  plane  one  discussed 
by  Benjamin.  Instability  either  by  divergence  or  by  flutter  is  pos- 
sible9but  due  to  the  presence  of  an  additional  parameter  which  measures 
the  "out-of-planeness"  of  the  system,  a  wider  variety  of  types  of  be¬ 
havior  is  possible. 

The  feature  of  the  spatial  system  chosen  is  that  one  can  experi¬ 
mentally  observe  the  various  modes  of  instability  by  simply  varying 
the  "out-of-planeness"  of  the  system.  3y  contrast,  in  Benjamin's  work, 
it  was  necessary  to  use  sets  of  pipes  of  different  mass  densities  to 
observe  different  types  of  instability. 


2.  The  System 

We  consider  here  a  generalization  of  Benjamin's  system,  allowing 
the  pipe  segments  to  oscillate  in  two  ditferent  planes.  The  system 
consists  of  two  straight  pipe  segments,  with  the  upper  segment  pinned 
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at  its  upper  end  0  in  such  a  manner  that  it  is  constrained  to  move  in 
a  single  vertical  plane,  the  x-y  plane  as  shown  in  Fig.  1.  The  lower 
segment  is  pinned  to  the  upper  one  at  P  such  that  an  angle  g  is 
formed  between  the  z-axis  and  the  normal  to  the  plane  of  motion  of  the 
lower  segment  (i.e.  the  pin  axis  e^  ).  Thus  for  3=0  the  motion  of 
the  two  segments  occurs  in  the  same  plane,  while  for  0  =  90°  the  planes 
of  motion  are  normal  to  one  another. 

An  incompressible  fluid  enters  the  upper  segment  at  0  and  is  dis¬ 
charged  at  the  free  end  of  the  lower  segment.  It  will  be  assumed  that 
the  rate  of  discharge  of  the  fluid  is  constant,  i.e.  that  the  fluid 
velocity  is  constant.  At  the  joints,  linear  restoring  springs  and  lin¬ 
ear  viscous  damping  are  introduced. 


3.  The  Equations  of  Motion 


Following  Benjamin  [5],  the  equations  of  motion  of  the  system  with 
two  degrees  of  freedom  under  consideration  may  be  derived  from  the  La- 
grangian  form  • 


d_  /3  L\  3L_  SD_ 
dt  \84i/  ‘  ^  aql 


3R 

P  VA  (V£t  +  R)  *  3=-  ; 


1=1,2 


where  q^  are  the  angles  shown  in  Fig.  1  and  q^  are  the  correspond¬ 
ing  angular  velocities  and 

R  =  Position  vector  of  free  end  (see  Fig.  1) 

=  Unit  vector  tangent  to  the  free  end 
V  =  Fluid  flowrate  relative  to  the  pipes 


Tl  + 


T  - 
a2 


Vl" 


V* 


-4- 


=  Kinetic  energy  of  the  pipes  alone 
T 2  =  Kinetic  energy  of  the  fluid  instantaneously  within  the  pipes 
=  Potential  energy  of  the  pipes  alone 
V2  =  Potential  energy  of  the  fluid  instantaneously  within  the  pipes 
A  =  Inside  cross-sectional  area  of  the  pipes,  assumed  the  same 
for  both  segments 
P  =  Mass  density  of  the  fluid 
D  =  Rayleigh's  Dissipation  function. 

In  the  following,  the  segments  are  assumed  to  be  axisymmetric  about  the 
center  line  of  the  flow  channel.  The  derivation  of  the  equations  and 
their  linearization  about  the  position  of  equilibrium  q^  a  q2  =  0 
is  straightforward.  The  potential  and  kinetic  energies  are  given  by 
V1  =  "  MA  g  a  cos  q1  -  0.5(1^  q^  +  K2  q22) 

-  Mg  g[(LA  +  b  cos  q2)  cos  q1  -  b  Cp  sin  qj^  sin  q 
v2  =  -  P  A  La  g(0.5  La  cos  qp 

-  p  A  Lg  g[(LA  +  0.5  Lg  cos  q2>  cos  q^  -  0.5  Lfi  Cp  sin  q^^  sin  q2] 

2Xl  -  +  Mg{[Sp  qt(LA  +  b  cos  q2)]2 

+  [b  q2  +  Cp  q^(b  +  LA  cos  q2>]2  +  [I*A  ^  Cp  sin  q2]2} 

+  ij^  Sp  cos  q2)2  +  (q2  +  Cp  q^2  I®,  +  Sp  sin  q^2!^ 


5  " 
\ 


k 


!• 
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2T2  -  pA|V2  La  +  La3  qL2/3  +  (q£  +  C0  q^2  Lfi3/3 

+  LB2(q2  +  cp  ^1>  La  ql  coS  q2J 
+  LgC^  La  Cp  cos  q2)2  +  Lb(La  ^  Sp)2 

+  IA  L02  cos  q2(<l1  SP)2  +  Lij3^1  cos  q2^3 

+  Lg  V2  +  2V  La  Lg  qt  CP  sin  q2  +  L0(LA  ^  CP  sin  q.,)2j- 

where 

Ik 

zz  =  Moment  of  inertia  of  upper  segment  with  respect  to  z  axis 
through  point  0  (see  Fig.  1) 

I®1  ,  I®2  ,  I33  ®  Principal  moments  of  inertia  for  the  mass  center 

of  the  lower  segment  with  respect  tc  e^  ,  je2  and  axes 

M.  ,  Md  =  Masses  of  upper  and  lower  segments 
A  B 

La  ,  Lg  «  Lengths  of  upper  and  lower  segments 
a  =  Distance  from  mass  center  of  upper  segment  to  point  0 

b  =  Distance  from  mass  center  of  lower  segment  to  point  P 

g  =  Acceleration  of  gravity 
CP  =  cos(p) 

Sp  =  sin(P"> 

K.^  ,  K2  =  Linear  spring  constants  at  upper  and  lower  joints  , 
respectively 

Rl  ,  R2  -  Linear  viscous  damping  constants  at  upper  and  lower 
joints, respectively 


The  equations  of  motion  are  then  obtained  as 


’i[4  •>  ifi  +  «*  b2  +  «B  la2  +  b  la 


+  p  A(La3  +  Lb?)/3  +  p  a  La  Lb(La  +  I,B)  ] 

+  +  Mg  b2  +  Kg  La  b  +  p  A  Lg3/3  +  p  A  La  Lb2/2 )] 

+  qx[p  A  V(La  +  Lb)2  +  Rl]  +  q2[Cp  p  A  V(2LA  Lg  +  Lg2)  ] 

+  q],[MA  ga  +  Mg  g(LA  +  b)+0.5pA  g(LA  +  Lg)2  + 

+  q2[CP(0.5  p  A  g  Lg2  +  Mj  g  b  +  D  A  V2  LA)3  =  0 


^[03(1^  +  Mg  b2  -f  Mg  b  La  +  p  A  Lb3/3  +  p  A  La  Lg2/2)] 
+  qjtlJj  +  Mg  b2  +  p  A  LR3/3] 

+  q^CB  p  A  V  Lg2]  +  q2[P  A  V  Lg2  +  R2] 


+  q1[C3(0.5  p  A  g  LgZ  +  Mg  g  b)] 

+  q2[0.5  p  A  g  Lg2  +  Mg  g  b  +  K2]  ■=  C 


In  the  first  part  of  the  paper  attention  is  restricted  to  the 
idealized  case  in  which  the  upper  and  lower  segments  are  straight  pipes 
of  uniform  cross-section  and  no  spring  or  damping  forces  act  at  the 
joints.  For  this  case: 


K1  -  K2  =  0 


*1  -  B2  -  0 


— ,  „  tjff^  -1 


**’?■,-'*»’»■  v*“  w  t 
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a 


b 


“  Mf  la2/3 

-  «B  LB2/12 

-  V2 

'  V2  , 


Following  Benjamin  [Si, we  assume  the  pipes  to  have  the  same  mass  density 
per  unit  length,  r;  .  Then ,introducing  the  following  non-dimensional 
parameters 


*\ 


a  =  la/lb 


V  =  3  p  A/(n  +  pA) 
u  -  Y  V/  Jl.S  g  Lb 

*k 

t  =/Ts  g/L  t 


(3) 


the  equations  of  motion  are  transformed  to  the  non-dimensional  form 


(<*  + 1)3 

0.5  C.p(2  +  3a) 

0.5  C3(2  +  3a) 

1 

ql" 

Lq2n 

- 

— 

L  J 

(Qf  +  l)2  U 

cp  u 

cp(l  +  2»)  U 

U 

V 

S- 

- 

— 

(a  +  l)2 

C3 

Cpj  (1  +  q  U^/y) 

1 

i 

ql 

_q2_ 

where  (  )'  denotes  d/dt 


(4) 
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4.  Stability  Analysis 

We  shall  examine  the  nature  of  the  solutions  of  the  form 

qi  =  Qi  eXt  ’  i=1’2  (5) 

which,  upon  substitution  into  the  equations  of  motion,  leads  to  the 
characteristic  equation 

/  ^  2 

A^  X  +  A^  X  +  A£  X  +  A^  X  +  Aq  =  0  (6) 

where 

A4  =  cy3  +  cy2(3  -  2.25  Cp2)  +  (1  +  3a)  Sp2 

A,  =  U[a3  +  q-2(4  -  3Cp2)  +  5o-  Sp2  +  2Sp2] 

A2  =  (q-  +  l)3  +  (cy  +  l)2  (1  +  U2) 

-  Cp2[(l  +  1.5  cy)  (2  +  q-U2/y)  +  (1  +  2  cy)U2] 

Al  =  2(cy  +  l)2  U  -  Cp2[2(cy  +  1)  U  +  cyU3/y] 

Aq  =  (cy  +  l)2  -  Cp2[1  +  cyU2/y] 


The  position  of  equilibrium  q^  =  q2  =  0  is  said  to  be  stable  provided 
all  the  characteristic  exponents  have  negative  real  parts.  The  well- 
known  Routh-Hurwitz  [3]  criterion,  applied  to  the  fourth  order  poly¬ 
nomial  at  hand,  assures  stability  if 

>  0  i=0 , 1 , . . . ,  4 

^  X  h  A^  A2  A^  ■  Aq  A^  -  A^  A^  >  0 


(8) 


r*W!Ni_'U' 
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Xf  any  of  the  above  conditions  are  violated,  the  system  will  be  un¬ 
stable. 

A  simplification  is  possible  since  both  and  A^  are  always 

positive.  For  small  U  ,  all  A^  are  positive.  Suppose  as  U  in¬ 
creases  A^  vanishes,  while  the  remaining  A^  are  positive.  Then 
X  must  be  negative,  and  hence  X  vanishes  before  A^  .  The  same 
applies  to  A^  .  Thus  the  six  requirements  for  stability  (8)  may  be 
replaced  by  just  two  conditions 


{ 


A0>° 


X  >  0 


(9) 


Loss  of  stability  can  thus  occur  in  two  different  ways.  If  Aq  =  0  , 
there  will  be  a  root  X  =  0  which  corresponds  to  an  adjacent  equilibrium 
position,  i.e.  to  static  Euler  buckling  (divergence).  By  contrast  when 
X  <  0  ,  the  characteristic  exponents  X  will  be  complex  with  positive 
real  parts,  and  stability  will  be  lost  by  oscillations  with  increasing 
amplitude  (flutter).  Written  out  we  have 

X  =  U6*[[-(d7y)(<v  +  1)V3  +  IV2  +  50-  +  2)  ]  Cp2 
+  CpT(<*/y)(a  +  l)2  (5a-2  +  10a-  +  4) 

+  (a/y )2(a  +  1)2(|  a2  v  3a  +  l)] 

+  Cp6(*-  (a/y)  (6a3  +  1  'x 2  1  9a  +  2) 

•  (o'/y)2^)  (18a3  +  33a2  +  20a  +  4)]  ]■ 

+  lAj^2[(a  +  1)^  +  (a  +  I)**] 
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•f  Cg2(a  +  l)4  C-2(6o'2  +  13a  +  6)  -  (a/Y)(3cy2  +  4a)] 

+  Cg4(a  +  l)2  [22a3  +  56a2  +  46a  +  12  +  (a-/y)(9a3  +  11a2  +  4a')] 

+  CP6  (a-  +  1)2[-2(2q-  +  1)(2  +  3a)]} 

+  U2^[a2(Q'  +  l)6]  +  Cg2[(a  +  l)4(2a2  -  2a3)] 

+  Cg4[-  (a  +  l)3(3a2  +  3a3)]}  (10) 

If  the  flowrate  U  is  very  small,  the  position  of  equilibrium  = 
q£  =  0  is  clearly  stable.  As  the  flowrate  is  slowly  increased,  Benja¬ 
min  [5]  has  shown  that  in  his  plane  case  (g  =  0)  instability  occurs 
first  by  buckling  if  y  >  0.5  and  by  flutter  if  y  <  0.5  .  It  is  found 
that  in  the  present  spatial  case  (g  ^  0)  the  value  Y  =  0.5  still 
separates  the  two  different  modes  of  loss  of  stability. 

The  Case  y  >  0.5 

It  is  proposed  to  show  that  if  y  >  0.5  ,  loss  of  stability  can 
occur  only  by  buckling,  as  the  flowrate  is  increased,  no  matter  what 
the  value  of  g  .  To  accomplish  this  it  must  be  shown  that  Aq  be¬ 
comes  zero  before  X  does,  as  the  flowrate  is  increased. 

The  flowrate  for  vanishing  Aq  ,  denoted  ,  is  obtained  from  (7) 
as 

uB2  =  (y/«)[(c*  +  i)2/ce2  -  i]  (id 

As  g  increases,  Ug  ir.  creases  and  at  g  =  90°,  no  finite  value  of  U 
can  cause  instability. 

On  substitution  of  Ug  into  the  expression  for  X  given  by  (10), 
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a  polynomial  in  powers  of  a  is  obtained  of  the  form 

9 

X  -  [VV  Cp2)]  )  ®n[Bnl  y.+  Bn2]  (12) 

n=l 

9 

where  the  are  polynomials  in  Cp  .  The  coefficient  of  a  is 

given  by 

V  -  0.5  Cp2 


and  thus  it  is  seen  that  a  necessary  condition  for  X  to  be  positive 
for  all  a  and  B  is  that  y  >  0.5.  Further,  it  is  readily  shown  that 
the  terms  are  all  non-negative  in  the  remaining  coefficients. 

Hence  the  value  of  X  for  a>0.5  must  be  larger  than  the  value  of 
X  for  a  «*  0.5.  Evaluating  the  latter  one  obtains 


[X] 


y=0.5 

tMJB 


[ub2/(o  CP2)]  jor8[2.5(l  +  sp2)  S32] 

+  a7  [(4/3  -  CP2)(16.5  -  18.75  Cp2  +  3Cp4)] 

+  O'6  [Sp2(56  -  67.5  CP2  +  18  CP4)] 

+  a-5  [Sp2(91  -  135.5  Cp2  +  47  Cp4)] 

+  o-4  [Sp2 (98  -  74.5  Cp2  +  1.5  Cp4)] 

+  a3  [sp2(70  -  72.5  Cp2  +  7.75  Cp4)] 

+  a2  [SP2(32  -  42.5  Cp2  +  10.5  Cp4)] 

+  [SP8]  }  (13) 


It  is  observed  that  the  coefficient  of  each  power  of  a  above  is  non¬ 
negative.  Hence  at  U  *»  Ug  ,  X  is  positive  for  y  >  1/2  . 
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2  2 

From  (10)  it  is  seen  that  X/U  is  a  quadratic  in  U  ,  i.e. 

X/U2  -  A  U4  +  B  U2  +  C 

and  also  that  C  is  always  positive.  For  small  U  ,  clearly  X  >  0 
for  all  y  .  Further,  for  0.5  <  y  <  3  it  is  found  that  X  >  0  in 
the  whole  interval  0  <  U  <  .  It  is  thus  concluded  that  if  y  >  0.5 

loss  of  stability  occurs  only  by  ouckling  as  the  flowrate  is  increased 
from  a  value  of  zero. 

It  might  be  of  interest  to  examine  whether,  at  flow  velocities 
higher  than  ,  a  transition  to  stability  or  flutter  could  occur. 

A  numerical  investigation  was  carried  out  for  the  following  range  of 
parameters 

0.1  <  a  <  10 

0.5  <  y  <  3 

0  "Z  U  <  10  uB 

For  numerous  values  of  a  ,  g  and  y  the  roots  of  (6)  were  computed 

as  U  increased  from  zero.  A  typical  root  locus  plot  is  shown  in  Fig.  2 
corresponding  to  g  =  0°  .  Since  any  complex  roots  must  occur  as  con¬ 
jugate  pairs,  only  the  upper  half  of  the  plane  is  shown.  When  U  =  0  , 
the  characteristic  roots  are  all  purely  imaginary,  and  as  U  is  in¬ 
creased,  the  roots  trace  out  the  curves  shown.  When  U  =1.56  a  pair 
of  roots  coalesces  on  the  negative  real  axis,  and  as  U  is  further  in¬ 
creased,  they  split,  one  moving  in  the  positive  and  one  in  the  negative 
real  directions.  Buckling  first  occurs  with  the  vanishing  of  one  root 
when  U  =  1.732  and  for  all  greater  U  ,  at  least  one  root  is  positive 
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real.  As  seen  in  the  figure,  all  roots  are  purely  real  after  U  =  3.12  , 
and  au  U  is  increased,  they  monotonically  increase  in  magnitude.  Thus 
divergence  is  the  only  mode  of  loss  of  stability  possible  for  y  >  0.5  t 
even  for  U  greater  than  Ug  . 

The  Case  v  <  0.5 

This  case  is  more  complex  in  that  either  buckling  or  flutter  may 
occur  first  as  the  flowrate  is  increased,  depending  on  the  value  of  3  • 
For  3  c  0  ,  Benjamin  [5]  has  shown  that  stability  will  be  lost 
first  by  flutter.  For  3  -  90°,'  on  the  other  hand,  * 

X  =  U4[2(o-  +  l)7  +  2(o-  +  l)6]  +  U2[q-2(q-  +  l)6] 

which  is  always  positive,  and  thus  is  positive  also  in  some  neighborhood 

of  $  =  90°  .  Since  Aq  can  vanish  for  all  3  except  8  °  90°,  it  is 

seen  that  near  3  =  90”  only  buckling  can  occur.  Thus  as  9  increases 

from  zero,  a  transition  in  the  mode  of  loss  of  stability  from  flutter  to 

* 

buckling  has  to  occur  at  some  value  of  the  angle,  3 

To  determine  this  transition  angle  3  ,  and  to  obtain  quantitative 
data  of  the  critical  value  of  U  for  the  whole  range  0  <  $<  9Cf,  a  param¬ 
etric  computer  study  was  carried  out.  A  typical  plot  of  U  as  a  func¬ 
tion  of  3  for  a  given  system  (a  «=  1.0  ,  V  =  0.25)  is  shown  in  Fig.  3. 
In  this  case  3  =26°.  In  region  I  there  are  two  pairs  of  complex  con¬ 
jugate  roots,  one  of  which  has  positive  real  parts.  Hence  in  this  region, 
loss  of  stability  occurs  by  flutter.  In  region  II  there  are  s  pair  of 
flutter  roots  and  two  real  roots,  one  of  which  is  positive.  Hence  in 
this  region  the  system  experiences  loss  of  stability  by  osci  .lation 
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with  increasing  amplitude  superimposed  upon  a  monotonic  motion  away  from 
the  equilibrium  position.  In  region  III  the  characteristic  roots  are 
all  real,  with  at  least  one  being  positive.  Hence  this  is  a  region  of 
divergence . 

A  typical  root  locus  plot  for  (3  =  0°  is  presented  in  Fig.  4.  As 

before  only  the  upper  half  of  the  plane  is  shown.  As  in  the  case  y  > 

0.5  ,  once  all  the  roots  become  real,  they  monotonically  increase  in 

magnitude.  Hence  there  are  no  new  instability  regions  above  the  buckling 

* 

region  in  Fig.  3.  For  all  3  >  3  >  the  typical  root  locus  is  the  same 

as  that  shown  in  Fig.  2  for  the  case  y  >  0.5  . 

* 

A  composite  plot  of  the  dependence  of  3  on  y  and  a  is  pre¬ 
sented  in  Fig.  5.  For  any  given  physical  configuration  (i.e.  specified 
values  of  a  ,  3  and  Y),one  can  determine  first  the  mode  of  loss  of 
stability.  Next,  the  corresponding  critical  f lowrate  can  be  calculated 
from  either  (10)  or  (11).  It  is  noted  that  as  y  — >  0.5  ,  3  — •>  0 
independent  of  a  ,  which  is  consistent  wi'th  the  ‘results  of  the  case 
Y  >  0.5  . 


5 .  Experiments 

An  overall  view  of  the  experimental  set-up  is  given  in  Fig.  6.  The 
tube  segments  were  of  thin  metal  (copper  and  brass)  in  commercially 
available  standard  gauges,  and  were  of  equal  lengths. 

The  joints  consisted  of  ball  bearings  and  light  brackets  made  of 
plexiglas.  The  pipes  were  connected  across  the  joints  by  short  segments 
of  latex  surgical  tubing  (Kent  Latex  Products,  0.25  inch  ID,  Thin  Wall). 
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Water  was  used  as  the  fluid,  and  its  velocity  was  measured  by  means  of 
a  precision  flowmeter  of  the  rotameter  type. 

The  additional  masses  of  the  hinges,  while  small,  were  not  negli¬ 
gible.  Further,  it  was  found  that  the  surgical  tubing  gave  rise  to 
small  restoring  and  damping  forces.  A  schematic  of  the  apparatus  is 
shown  in  Fig.  7.  The  quantities  needed  in  the  governing  equations  can 
be  identified  as 

M,  «  ra  + 

A 

Mg  =  m  +  i 

LA  =  x  + 

h°i  + 1 

MA  a  =  (m 

Mg  b  =  mC 

A  m.C,^ 

Izz  ”  TT 

*U  +  V 

where 

m  c  mass  of  each  metal  tube 
l  =  length  of  each  metal  tube 
m^  B  mass  of  each  hinge 

Ir  “  In  of  each  hinge  with  respect  to  its  centroid. 


“h 

2c 

c 

+  2mg)(0.5j&  +  c) 

0.5X  +  c)  +  nijj  d 

+  m(0.5i  +  c)2  +  2^  +  mg  d2  +  1^(4  +  2c  -  d)2 
2 

b2  =  ||-  +  m(0.5X  +  c)2  +  In  +  d2 
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Xt  was  suspected  that  the  restoring  and  damping  forces  would  be 
sensitive  to  the  water  pressure  in  the  surgical  tufing  at  the  joints. 

The  values  of  the  restoring  and  damping  constants  were  obtained  by 
conducting  a  separate  test  with  just  the  upper  tube  and  upper  hinge. 

The  tube  was  closed  at  its  lower  end  and  pressurized  with  air.  It 
was  then  photographed  with  a  motion  picture  camera  as  it  oscillated 
freely  and  damped  out.  This  was  repeated  for  several  different  pres¬ 
sures.  From  the  observations  the  restoring  and  damping  constants 
could  be  obtained  in  a  standard  fashion. 

To  utilize  this  information  it  was  necessary  to  know  the  fluid 
pressure  in  the  surgical  tubing  as  a  function  of  flowrate.  This  was 
obtained  by  assembling  the  pipes  and  hinges  in  the  test  configuration 
(as  in  Fig.  6)  and  replacing  one  joint  at  a  time  by  a  short  tube 
tapped  to  place  a  pressure  gauge.  Then  water  was  forced  through  and 
pressure  readings  could  be  taken.  This  procedure  neglects  the  effects 
of  expansion  of  the  surgical  tubing  and  of  motion  of  the  joints  with 
increase  in  pressure.  For  small  motions,  however,  the  latter  should 
have  little  effect,  and  the  final  tests  were  discontinued  when  the  ex¬ 
pansion  of  surgical  tubing  became  excessive. 

The  results  of  the  above  preliminary  tests  indicated  that  the  lower 
spring  constant,  ,  was  nearly  independent  of  the  flowrate  over  the 
range  involved,  while  the  upper  spring  constant,  ,  was  fairly  sen¬ 
sitive  to  the  flowrate,  as  shown  in  Fig.  8.  This  could  have  been  ex¬ 
pected,  since  most  of  the  pressure  drop  involved  would  take  place  at  the 
upper  joint  where  the  flow  channel  was  constricted  down  to  the  inside 
diameter  of  the  pipe.  It  was  also  found  that  the  damping  coeffic Lents 


were  nearly  independent  of  pressure,  and  thus  the  damping  is  primarily 
due  to  the  friction  in  the  ball  bearings. 

The  actual  experiments  to  determine  the  instability  mode  behavior 
were  performed  on  two  different  pairs  of  tubes,  one  with  y  =  0.328  and 
one  with  Y  =  1.09  .  The  dimensions  of  the  two  sets  of  pipes  and  of 
the  hinges  are  shown  in  the  table  below.  The  same  hinges  were  used  with 
each  set  of  tubes. 


Set  I  y  =  0.328 

Set  II  y  =  1.09 

Type 

Copper  l/4"0Dx.035  wall 

Brass  5/16"ODx  .015  wall 

m  (slugs) 

0.000198 

0.000123 

m^slugs) 

0.000059 

0.000059 

i  (inch) 

12.0 

12.0 

A  (eq.inch) 

0.0254 

0.0625 

K2  (in- lb) 

0.206 

0.206 

R.  ,R„(in-lb-sec) 

1  4 

0.030 

0.030 

c  (inch) 

0.750 

0.750 

d  (inch) 

-  ,  2 

1 . 295 

1.295 

IR  (in  -slug) 

0.000135 

i 

0.000135 

With  these  values,  the  coefficients  in  the  governing  equations  were 
evaluated  and  the  crifLcal  flowrates  obtained  as  described  in  the  first 
part  of  the  paper.  The  results  are  shown  in  Figs.  9  and  10.  The  solid 
curves  denote  the  theoretical  predictions. 

The  tests  were  performed  quasi-statically .  Thus  the  flowrate  was 
increased  in  small  increments,  and  at  each  step  the  two  pipes  were 
manually  displaced  two  degrees  from  the  vertical  and  released  with  no 


initial  velocity.  This  was  repeated  until  a  flowrate  was  reached  at 
which  the  resulting  motion  did  not  subside.  The  onset  of  flutter  was 
marked  by  steady  periodic  oscillations.  Upon  slightly  increasing  the 
flowrate,  there  resulted  oscillations  with  increasing  amplitude. 

In  certain  cases,  the  onset  of  buckling  was  somewhat  difficult  to 
determine  precisely.  This  was  due  to  unavoidable  small  eccentricities 
in  the  system  which  caused  the  position  of  equilibrium  to  change 
slightly  with  an  increase  in  flowrate.  As  a  result, the  onset  of  buck¬ 
ling  was  evidenced  by  large  deflections  corresponding  to  a  finite  (but 
narrow)  range  of  flowrates.  In  the  majority  of  cases,  however,  as  the 
flowrate  was  increased  in  small  steps,  the  transition  to  the  buckled 
state  occurred  quite  suddenly.  The  associated  deflections  were  so 
large  as  to  require  steps  to  prevent  the  system  from  being  damaged. 

Thus  in  most  cases  an  unambiguous  buckling  load  could  be  determined. 

The  above  procedure  was  repeated  for  each  increasing  value  of  6  . 

The  range  of  the  data  taken  was  limited  by  the  expansion  of  the  surgical 
tubing  at  large  flowrates. 

The  tubing  gave  rise  to  another  problem.  As  it  is  manufactured 
and  packed,  it  has  a  slight  permanent  curvature.  When  the  pipes  and 
joints  were  assembled,  the  surgical  tubing  w&j  arranged  such  that  the 
plane  of  motion  of  the  joint  was  normal  to  Kh-i  plane  of  curvature  of 
the  tubing.  Without  this  precaution  the  1  wo  metal  tubes  could  not  be 
aligned  even  in  the  absence  of  flow.  Thus,  6  was  changed,  it 

was  also  necessary  to  rearrange  the  surgical  ti.bi.  ?.  And  in  doing  this 
great  care  had  to  be  taken  to  ascertain  that  the  tubing  was  not 
stretched.  It  was  found  early  in  the  experimentation  that  even  a  small 
amount  of  stretch  .np,  resulted  in  considerable  changes  in  the  buckling  load 


The  surgical  cubing  also  gave  rise  to  a  notable  non- linear  ef¬ 
fect.  When  stability  was  lost  by  flutter  it  was  found  that  an  in¬ 
crease  of  the  initial  displacement  resulted  in  a  lower  critical  flow- 
rate.  This  implies  the  existence  of  either  non-linear  restoring  or 
damping  forces  of  the  "softening"  type  and, in  fact,  it  was  found  from 
the  dynamic  amplitude- time  data  (used  to  determine  the  spring  con¬ 
stants)  that  non-linear  damping  was  the  predominant  factor.  In  the 
case  of  such  "softening"  type  forces  one  cannot  invoke  the  theorem  due 
to  Lyapunov  as  in  [3],  which  asserts  that  the  instability  threshold 
for  a  non-linear  system  is  the  same  as  that  for  the  corresponding 
linearized  system.  Work  is  underway  to  evaluate  this  effect  in  a 
quantitative  way,  but  until  then  this  does  serve  to  emphasize  the  need 
for  employing  uniform,  small  initial  conditions  in  tests  of  this  type. 

The  results  of  the  tests  are  summarized  in  Figs.  9  and  10.  For 

ic 

y  =  0.328,  the  transition  angle  8  at  which  the  initial  instability 
changed  from  flutter  to  buckling  was  predicted  analytically  to  lie 
between  21°  and  22°.  In  the  test,  the  pin  angle  8  was  increased 
in  5 6  increments  and  as  seen  in  Fig.  10  the  changeover  occurred  be¬ 
tween  8  =  20°  and  3  =  25°.  No  closer  determination  was  attempted, 
since  the  inaccuracies  in  the  values  of  the  spring  and  damping  con¬ 
stants  rendered  any  further  refinement  meaningless.  The  transition 
angle  provides  a  measure  of  the  correlation  between  the  experimental 
system  and  the  idealized  one  considered  in  the  first  portion  of  the 
paper.  For  the  latter,  with  y  =  0.328  ,  Fig.  5  predicts  a  transition 
angle  of  20°  and  this  agrees  well  with  the  transition  angle  for  the 

*j(f 

experimental  system,  which  was  predicted  as  21°  <  g  <  22°  and 
measured  as  20°  <3  <  25°  . 
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In  each  of  the  cases,  the  measured  uuckling  loads  were  lower  than 
those  predicted,  as  would  be  expected  due  to  the  unavoidable  presence 
of  small  eccentricities.  The  agreement  in  each  case  tends  to  fall  off 
as  p  is  increased.  For  the  flutter  case  (v  =  0.328)  the  discre¬ 
pancy  varied  from  4%  at  8  =  0°  to  10%  at  3  =  40°.  T\  che  buckling 
case  (v  =  1.09)  the  discrepancy  varied  from  8%  at  b  =  0°  to  18%  at 
8  =  40° .  The  increasing  discrepancy  as  8  is  increased  is  probably 
attributable  to  the  normal  forces  applied  to  the  bearing  which  are 
present  in  increasing  magnitude  as  8  is  increased.  This  would  tend 
to  accentuate  any  eccentricities  in  the  joint  and  further  cause  beam- 
type  bending  in  the  tubes  themselves,  giving  rise  to  further  misalign¬ 
ment.  In  view  of  the  unavoidable  imperfections  present  in  any  mechani¬ 
cal  system  used  to  study  buckling  behavior,  the  agreement  between  the 
theory  and  the  experiment  obtained  herein  is  considered  to  be  acceptable. 


5.  Concluding  Remarks 


considering  ?.  special  system  of  articulated  pipes,  it  has  been 
shown  that  cne  various  cypes  of  loss  of  stability,  characteristic  of 
idealised  non- conservative  systems,  can  in  fact  be  obtained  in  a  phys¬ 
ically  realizable  model.  It  has  also  been  verified  that  the  assumption 
of  one- dimensional  flow  and  the  use  of  linearized  equations  of  motion 
are  adequate  to  determine  the  boundaries  of  the  various  regions  of 
instability. 

However,  an  important  aspect  of  such  systems  has  been  observed, 
which  must  undoubtedly  be  considered  in  applying  the  above-mentioned 
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assumptions  to  stability  problems  in  pipes  conveying  fluid.  As  found 
in  the  experimentation,  the  softening- type  non-linear  behavior  of 
the  flexible  joints  will  significantly  change  the  stability  charac¬ 
teristics  of  the  system.  Thus  a  careful  characterisation  of  the  non¬ 
linear  behavior  of  coupling  tubes  is  a  necessary  pie  recorisite  to  the 
practical  stability  analysis  of  such  devices  as  fluid- dr 1  ''er>  control 
systems. 
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